Phase diagram of Bose-Fermi mixtures in one-dimensional optical lattices 
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The ground state phase diagram of the one-dimensional Bose-Fermi Hubbard model is studied in 
the canonical ensemble using a quantum Monte Carlo method. We focus on the case where both 
species have half filling in order to maximize the pairing correlations between the bosons and the 
fermions. In case of equal hopping we distinguish between phase separation, a Luttinger liquid phase 
and a phase characterized by strong singlet pairing between the species. True long-range density 
waves exist with unequal hopping amplitudes. 
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Ultracold atomic gases loaded in an optical lattice 
offer unique possibilities to simulate many fundamental 
models of solid state physics. The system is extremely 
clean and there is unprecedented experimental control 
over the lattice period and the interaction strength. 
The prediction Q and the successful observation of the 
bosonic superfluid-Mott transition by Greiner et al. 
was the key breaktrough and led to a revival in interest 
in bosonic physics. Experiments on fermions in different 
hyperfine states are under way Q and can potentially be 
employed for studies of high-Tc superconductivity and 
Mott-insulating phases. Ultracold atoms can also be ap- 
plied to create a completely different and unique state of 
matter: bosonic and fermionic atoms can simultaneously 
be trapped in order to let them interact in a controlled 
way. Nowadays, several stable Bose-Fermi mixtures 
have been created , e.g., 
87 Rb - 40 K H and 87 Rb 
an optical lattice. These experimental advances have 
challenged theoreticians: the prop erties of mixtures 
have been studied in Refs. 

using a variety of methods. Most remarkably, it was 
shown that exotic pairing between the bosons and the 
fermions might exist. However, the ground state phase 
diagram remained largely unknown, complicating any 
study to the influence of the experimentally unavoidable 
drawbacks, namely finite temperature and parabolic 
confinement. Using quantum Monte Carlo (QMC) 
simulations, we provide the complete, exact ground 
state phase diagram in one dimension at the parameters 
of greatest interest: double half filling with equal 
hopping tp = is • Afterwards, we discuss the case of 
strong anisotropic hopping tp = 4ta. Recently, the 
grand-canonical QMC simulations of Ref. focused 
on the issue of the stability of the mixed phase. 



We assume that a mixture of bosons and fermions is 
contained in a square well and loaded into an optical lat- 
tice. The temperature is low enough such that quantum 
degeneracy is achieved. The system is then described by 



a lowest-band Bose-Fermi Hubbard chain, 
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where and Cj are the bosonic and fermionic annihi- 
lation operators on site i, respectively. The intersite 
spacing a is set to unity, i.e., a = 1. The operators 
n i = Kbi and = c\ci denote the number operators 
on site i for bosons and fermions, respectively. Bosons 
(fermions) can hop from site i to a nearest neighbor site 
j = i±l with tunneling amplitude t& (tp). Furthermore, 
a large occupation of bosons on a single site is suppressed 
by the on-site repulsion term Ubb- Bosons and fermions 
can mutually repel or attract each other on every site 
depending on the sign of Ubf- In order to enhance 
the pairing correlations in the system, we assume that 
bosons and fermions both have a density corresponding 
to double half filling, N B = N F = L/2. We work in 
units tB = f, and simulate at inverse temperature 
j3 = 2L. In a homogeneous system at half fermionic 
filling, the particle-hole transformation, q — ► (— l) 4 c| 
changes the sign of the boson-fermion interaction in 
eq.JTJ, Ubf — > — ^bf- We choose Ubf > without loss 
of generality. 

One-dimensional phases can be characterized by the 
slowest decaying modes of the system. The bosonic 
Green function Gb(£ — j) = (b\bj), fermionic Green func- 
and composite pair Green func- 
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tion Gf(« — j) = 
tion Gbf(* — j) 
j) = (b\cjbjCi) (Ubf > 0) allow for such a classi- 
fication. In addition, the bosonic superffuid density 
Ps,b = WqL/(2(3) 17], the fermionic stiffness p s ,F = 
WpL/(2f3), the paired superffuid density ppsf = (Wb + 
Wf) 2 L/(2(3) [18| and the counter-rotating superfluid 
density p S CF = {W B - W F ) 2 L/(2P) [Tg| probe for soft 
modes in the system. The winding numbers for the 
bosons and the fermions are denoted by Wb and Wp, 
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respectively. The Green functions can either decay expo- 
nentially or by power law according to (for a continuum 
coordinate x 3> 1): 

G B (x)~(d(x\L))- K *, (2) 
G F (x) ~ sin{nx/2)d{x\L)- KF , (3) 
Gbf{x) ~ sm(Trx/2)d(x\L)- K * F 1 (4) 

where the cord function d(x\L) = | sin(7ra:/L)| takes care 
of the periodic boundaries |2jj . The relevant correlators 
for density-density correlations are Ccdw(Sdw) — 
{{m ± mi)(n ± m.j)) ~ ((n, ± m l ))((n J ± rrij)) ~ 
cos(7ra:)d(2;|L) _ cnw ( SDW ) , where the upper signs refer 
to the charge-density wave correlator (CDW) and the 
lower signs to the spin-density wave (SDW) correlator. 

We study the Bose-Fermi Hubbard model using a 
quantum Monte Carlo method based on worm-type 
updates |2l| , usin g a generalization of the canonical 
methods of Refs. [22, |23|- In a molecular superfluid 
phase bosons and fermions are paired into molecules 
and thus have to wind together. The standard single- 
particle worm algorithm changes the winding of the 
species independently. In order to create a molecular 
winding number, it needs to go through a state with 
different bosonic and fermionic winding numbers, which 
is exponentially suppressed. A key ingredient is the use 
of a two-particle, boson- fermion worm [24|, which can 
let the bosons and fermions move together. Technical 
details will be presented elsewhere. 

The main result of this study is the phase diagram 
of the Hamiltonian at double half filling and with 
parameters = t-p = t shown in Fig. ^ We distin- 
guish three phases: (i) Luttinger liquid phase (LL), 
(ii) (pseudo) spin-density wave phase (SDW) and (iii) 
phase separation (PS). The Green functions in the 
various phases are shown in Fig. [2 We now make a tour 
clockwise around the phase diagram, starting at {7bf = 0. 

When there is no interaction between the bosons and 
the fermions, Ubf — 0, we have a system consisting of 
free fermions and a Luttinger liquid of bosons. For small 
Ubf the fermions also behave as a Luttinger liquid, 
which weakly interacts with the bosonic liquid. The 
system thus consists of two weakly interacting Luttinger 
liquids (LL phase). The bosonic Green function is 
always the slowest decaying mode. For small values of 
Ubf, the fermionic Green function decays slower than 
the composite paired Green function (LLa), while there 
is a cross-over at larger U B f when the fermionic Green 
function decays faster (LLb). 

In the limit of infinite Ubb , the bosons behave as 
hard-core or Tonks bosons |2a , l27|b We can associate 
a pseudo-spin 'up' with the bosons and a pseudo-spin 
'down' with the fermions. The charge sector becomes 
SU(2) symmetric and the system can be mapped onto a 
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FIG. 1: Phase diagram of the Bose-Fermi Hubbard model, 
Eq.(0. The different phases are (i) Phase separation ('PS'), 
(ii) Luttinger liquid phase ('LLa' and 'LLb'), and (iii) a 
(pseudo) spin-density wave ('SDW'). The phase transition 
(solid line) to the PS phase is first order, while the tran- 
sition from LL to SDW belongs to the Kosterlitz-Thouless 
universality class |2S||. The dashed line indicates a cross-over 
when the composite pair Green function decays slower than 
the fermionic Green function. 
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FIG. 2: Bosonic('+', full line), fermionic('o', dashed line) and 
composite pair Green function('x', dotted line) in the phases 
LLa, LLb, SDW of Fig. Only the exponential decay or 
power-law behavior in absolute value of the correlation func- 
tions has been shown. The PS phase has been identified by a 
peak in the bosonic structure factor at momentum k — 2n/L, 
shown in the fourth plot. 



Heisenberg chain with effective spin-exchange amplitude 
J ~ 4£ 2 /[/bf for repulsive, small Ubf- The Heisenberg 
chain has only one phase, a spin-density wave (SDW) 
with a gapless spin sector. There is a charge gap which 
is exponentially small, A c ~ ^/UB~Fe.~ v ^ UsF , with v 
the velocity. In Fig. ^ we see that the SDW phase 
extends to finite values of Ubb- The phase is identified 
by exponential decay of the bosonic and the fermionic 
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Green functions, shown in Fig. EI The bosonic superfluid 
density and fermionic stiffness are small but non-zero, 
the pair superfluid density ppsf is zero for repulsive Ubf 
while pscf = 2(p s .B + Ps,f)- This is a reflection of the 
gapped charge and the gapless spin sector we identified 
in the limit Ubb — > +00. The transition line between the 
LL and SDW phases converges logarithmically slowly 
to zero, as a consequence of the exponentially small 
charge gap. In case of attractive Ubf, the SDW phase 
in Fig. ^ is replaced by a singlet paired phase, where the 
power-law decay of the composite pair Green function 
is slower than the decay of the charge-density wave 
correlator. The Pscf is zero and ppsf = 2(p s ,e + Ps,f)- 
When Ubb — > 00 the system becomes isotropic and the 
composite pair Green function decays as fast as the 
charge-density wave correlator. 

When both Ubb and Ubf are large, the 
model can be mapped onto an XXZ-Hamiltonian 
H XXZ = E« J(of of+i + ^f+i) + J z °l°l+i with 

J = -|^e and J z = - jf^ HI. Usin g the exact 

solution of the XXZ-chain, the first order transition 
towards phase separation (i.e., the ferromagnet) occurs 
at J = J z or Ubf /Ubb = 2 for ie = £f = 1. This is the 
slope of the transition line between the SDW and PS 
phases in Fig. ^ 

At infinite repulsion between the bosons and the 
fcrmions, the bosons and the fermions phase separate 
(PS) and form non-overlapping domains. When Ubb = 
and when Ubf is very large, the ideal Bose and Fermi 
gases are in first order unstable to phase separation 
with hard domain walls. The energy is minimized when 
the bosons occupy a number of y — 7r(L/2) 1 / 3 sites. 
This instability persists whenever Usp/t > 0(L 2 / 3 ) in 
first order when Ubf — ► 0. At larger Ubb, the boson 
repulsion exerts a pressure such that the region occupied 
by the bosons will grow. In Fig. |2 we see a clear 
peak at momentum k — in the bosonic structure 
factor. In the simulations we detect phase separation 
by a diverging peak in the structure factor at small 
momenta, accompanied by an exponential decay of 
all three Green functions and zero winding numbers, 
Ps,b = Ps,f = Psgf = Ppsf = 0. The phase transition to 
the phase separated phase is strongly first order. 

The physics of the different phases can also be 
understood when looking at the momentum profiles, 
which are experimentally accessible. In the LLa region, 
the fermionic momentum profile is very close to the 
non-interacting one. In the LLb, the fermionic momen- 
tum profile begins to deviate slightly, but the Fermi 
momentum is still well defined. In the SDW phase 
however, the fermions are strongly interacting while we 
see a Fermi surface for the paired composite particles. 
This is visualized in Fig. |5J 

Difference in masses between the bosons and the 
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FIG. 3: Fermionic momentum profiles in the LLa, LLb and 
SDW phases for the same system parameters as in Fig[5] The 
error bars are smaller than the point sizes. In the LL phase 
the Fermi momentum remains well defined, while the SDW 
phase destroys the Fermi surface of the bare fermions. The 
composite paired particles have a well defined Fermi momen- 
tum. 
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FIG. 4: Phase diagram of the Bose-Fermi Hubbard model, 
Eq.(0, at strong unequal hopping amplitudes. The differ- 
ent phases are (i) Phase separation (PS'), (ii) Luttinger liq- 
uid phase ('LL'), (iii) a gapped 'Neel (Ising)' state with true 
long-range order and (iv) a gapless (pseudo) spin-density wave 
('SDW') without long-range order. 



fermions will in general lead to strong unequal hopping 
amplitudes |2!j . We therefore examine the system with 
parameters tp = 4, £b — 1 at double half filling. At 
large Ubb and Ubf the XXZ Hamiltonian gives the 
asymptotic slope Ubf /Ubb — 4.5 between a gapless 
phase and an anti- ferromagnetic phase (Neel phase), and 
for the transition to the ferromagnet (phase separation) 
the slope is Ubf /Ubb — 12.5. The full phase diagram is 
shown in Fig.0] The Neel (Ising) state is gapped and has 
true long-range order. The density-density part of the 
Ising correlator converges to a constant, as is shown in 
Fig. [5j Note that the phase transitions happen at much 
larger values of Ubb and Ubf compared to the case of 
equal hopping. In the gapless SDW phase, the composite 
pair boson-hole Green function has the slowest decay, 
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FIG. 5: Decay of the composite pair (boson-hole) Green 
function (Gbf, '+' mark) and SDW order correlation 
function(CsDW, '*' mark) in the Neel state. The data points 
are interpolated by power-law fits, see Eq.0. 

while the pair superfluid density is zero, ppsf = 0. This 
is similar to the SDW state in Fig. ^ In case of Tonks 
bosons, Ubb — * oo, the transition from the LL phase to 
the Neel state belongs to the Kosterlitz-Thouless uni- 
versality class |3fJ. The transition from PS to the Neel 
phase cannot occur directly, but a gapless phase exists 
between PS and the Neel phase. The LL phase again 
has several crossovers, which we do not discuss here. 
In the limit Ubb — 0, Ubf — > the same perturbation 
argument as in the equal hopping case would again lead 



to phase separation in the thermodynamic limit. From 
the present numerics on small system sizes it is difficult 
to confirm this, although it cannot be excluded either 
and we see an inflection point in the transition line 
between PS and LL at low values of Ubf in Fig. 0] 

The phases presented in this paper can experimentally 
be identified as follows. The momentum profiles of the 
bosons, fermions and composite particles distinguish the 
weakly interacting phases from phases characterized by 
strong pairing, as shown in Fig. Noise correlations in 
time-of-fiight images can be used to find the gapless SDW 
and gapped Ising/Neel phases 0,H2- In principle, RF 
spectroscopy can be used to measure the energy gaps |33| . 

In conclusion, we have studied the ground state phase 
diagram of the Bose-Fermi Hubbard model at double 
half filling. Interactions can lead to such phases as 
a Luttinger liquid, a (pseudo) spin-density wave, a 
Neel(Ising) state and even phase separation can occur. 
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